Ozone Monitoring Instrument Observations of Interannual Increases in SO2 Emissions from Indian Coal-fired Power Plants During 2005-2012

Due to the rapid growth of electricity demand and the absence of regulations, sulfur dioxide (SO2) emissions from coal-fired power plants in India have increased notably in the past decade. In this study, we present the first interannual comparison of SO2 emissions and the satellite SO2 observations from the Ozone Monitoring Instrument (OMI) for Indian coal-fired power plants during the OMI era of 2005-2012. A detailed unit-based inventory is developed for the Indian coal-fired power sector, and results show that its SO2 emissions increased dramatically by 71 percent during 2005-2012. Using the oversampling technique, yearly high-resolution OMI maps for the whole domain of India are created, and they reveal a continuous increase in SO2 columns over India. Power plant regions with annual SO2 emissions greater than 50 Gg year-1 produce statistically significant OMI signals, and a high correlation (R equals 0.93) is found between SO2 emissions and OMI-observed SO2 burdens. Contrary to the decreasing trend of national mean SO2 concentrations reported by the Indian Government, both the total OMI-observed SO2 and average SO2 concentrations in coal-fired power plant regions increased by greater than 60 percent during 2005-2012, implying the air quality monitoring network needs to be optimized to reflect the true SO2 situation in India

A TROPOMI- and GLM-Based Estimate of NOx Production by Lightning over the U.S.

Lightning produces NO because the extreme temperatures (>20000 K) in lightning channels dissociate molecular O2 and molecular N2, which then combine to form NOx which quickly reacts with O3 to form NO2. Lightning is responsible for 10-15% of NOx emissions globally. This is 2 8 Tg N a-1 [Schumann and Huntrieser, 2007] or 100 to 400 mol per flash. Much of the uncertainty stems from limited knowledge of lightning NOx production per flash (LNOx PE) or per unit flash length. Most LNOx is injected into mid- and upper-troposphere where away from deep convection its lifetime is longer relative to lower troposphere NOx. NOx in this region enhances the concentrations of upper tropospheric NOy, OH, and O3 and contributes to positive radiative forcing by O3 and negative forcing by CH4. We have previously used OMI NO2 to obtain estimates of LNOx production per flash over the Gulf of Mexico (Pickering et al., 2016, JGR), in convective events during NASAs TC4 field program (Bucsela et al., 2010, JGR), and over broad regions of the tropics (Allen et al., 2019, JGR) and midlatitudes (Bucsela et al., 2019, JGR). In the latter studies, we obtained PE values of 170 100 mol flash and 180 100 mol flash, respectively

Effect of Tropospheric Aerosols in Satellite-Based Trace Gas Retrieval

Scattering and absorption by tropospheric aerosol particles have an effect on the airmass factor (AMF), a fundamental quantity for trace gas concentration retrieval by inversion of satellite measurements. The interference effect depends on the aerosols micro-physical and optical properties as well as the relative distribution of the tropospheric trace gas and aerosol load. Quantifying the aerosol impact on trace gas retrieval requires a sensitivity study using radiative transfer calculations. In this presentation we will describe a recently initiated effort to characterize the aerosol-related error in trace gas retrievals when the presence of aerosol particles is not accounted in the inversion procedure. A general description of the project will be presented and preliminary results on aerosol effects on SO2 retrieved concentrations will be discussed

New-generation NASA Aura Ozone Monitoring Instrument (OMI) volcanic SO2 dataset: Algorithm description, initial results, and continuation with the Suomi-NPP Ozone Mapping and Profiler Suite (OMPS)

Since the fall of 2004, the Ozone Monitoring Instrument (OMI) has been providing global monitoring of volcanic SO2 emissions, helping to understand their climate impacts and to mitigate aviation hazards. Here we introduce a new-generation OMI volcanic SO2 dataset based on a principal component analysis (PCA) retrieval technique. To reduce retrieval noise and artifacts as seen in the current operational linear fit (LF) algorithm, the new algorithm, OMSO2VOLCANO, uses characteristic features extracted directly from OMI radiances in the spectral fitting, thereby helping to minimize interferences from various geophysical processes (e.g., O3 absorption) and measurement details (e.g., wavelength shift). To solve the problem of low bias for large SO2 total columns in the LF product, the OMSO2VOLCANO algorithm employs a table lookup approach to estimate SO2 Jacobians (i.e., the instrument sensitivity to a perturbation in the SO2 column amount) and iteratively adjusts the spectral fitting window to exclude shorter wavelengths where the SO2 absorption signals are saturated. To first order, the effects of clouds and aerosols are accounted for using a simple Lambertian equivalent reflectivity approach. As with the LF algorithm, OMSO2VOLCANO provides total column retrievals based on a set of predefined SO2 profiles from the lower troposphere to the lower stratosphere, including a new profile peaked at 13 km for plumes in the upper troposphere. Examples given in this study indicate that the new dataset shows significant improvement over the LF product, with at least 50% reduction in retrieval noise over the remote Pacific. For large eruptions such as Kasatochi in 2008 (∼1700 kt total SO2/ and Sierra Negra in 2005 (\u3e 1100DU maximum SO2/, OMSO2VOLCANO generally agrees well with other algorithms that also utilize the full spectral content of satellite measurements, while the LF algorithm tends to underestimate SO2. We also demonstrate that, despite the coarser spatial and spectral resolution of the Suomi National Polar-orbiting Partnership (Suomi-NPP) Ozone Mapping and Profiler Suite (OMPS) instrument, application of the new PCA algorithm to OMPS data produces highly consistent retrievals between OMI and OMPS. The new PCA algorithm is therefore capable of continuing the volcanic SO2 data record well into the future using current and future hyperspectral UV satellite instruments

Airborne MAX-DOAS Measurements Over California: Testing the NASA OMI Tropospheric NO2 Product

Airborne Multi-AXis Differential Optical Absorption Spectroscopy (AMAX-DOAS) measurements of NO2 tropospheric vertical columns were performed over California for two months in summer 2010. The observations are compared to the NASA Ozone Monitoring Instrument (OMI) tropospheric vertical columns (data product v2.1) in two ways: (1) Median data were compared for the whole time period for selected boxes, and the agreement was found to be fair (R = 0.97, slope = 1.4 +/- 0.1, N= 10). (2) A comparison was performed on the mean of coincident AMAX-DOAS measurements within the area of the corresponding OMI pixels with the tropospheric NASA OMI NO2 assigned to that pixel. The effects of different data filters were assessed. Excellent agreement and a strong correlation (R = 0.85, slope = 1.05 +/- 0.09, N= 56) was found for (2) when the data were filtered to eliminate large pixels near the edge of the OMI orbit, the cloud radiance fraction was2 km, and a representative sample of the footprint was taken by the AMAX-DOAS instrument. The AMAX-DOAS and OMI data sets both show a reduction of NO2 tropospheric columns on weekends by 38 +/- 24% and 33 +/- 11%, respectively. The assumptions in the tropospheric satellite air mass factor simulations were tested using independent measurements of surface albedo, aerosol extinction, and NO2 profiles for Los Angeles for July 2010 indicating an uncertainty of 12%

Version 2 Ozone Monitoring Instrument SO2 product (OMSO2 V2): New anthropogenic SO2 vertical column density dataset

The Ozone Monitoring Instrument (OMI) has been providing global observations of SO2 pollution since 2004. Here we introduce the new anthropogenic SO2 vertical column density (VCD) dataset in the version 2 OMI SO2 product (OMSO2 V2). As with the previous version (OMSO2 V1.3), the new dataset is generated with an algorithm based on principal component analysis of OMI radiances but features several updates. The most important among those is the use of expanded lookup tables and model a priori profiles to estimate SO2 Jacobians for individual OMI pixels, in order to better characterize pixel-to-pixel variations in SO2 sensitivity including over snow and ice. Additionally, new data screening and spectral fitting schemes have been implemented to improve the quality of the spectral fit. As compared with the planetary boundary layer SO2 dataset in OMSO2 V1.3, the new dataset has substantially better data quality, especially over areas that are relatively clean or affected by the South Atlantic Anomaly. The updated retrievals over snow/ice yield more realistic seasonal changes in SO2 at high latitudes and offer enhanced sensitivity to sources during wintertime. An error analysis has been conducted to assess uncertainties in SO2 VCDs from both the spectral fit and Jacobian calculations. The uncertainties from spectral fitting are reflected in SO2 slant column densities (SCDs) and largely depend on the signal-to-noise ratio of the measured radiances, as implied by the generally smaller SCD uncertainties over clouds or for smaller solar zenith angles. The SCD uncertainties for individual pixels are estimated to be 0.15-0.3DU (Dobson units) between 40 S and 40 N and to be 0.2-0.5DU at higher latitudes. The uncertainties from the Jacobians are approximately 50 %-100% over polluted areas and are primarily attributed to errors in SO2 a priori profiles and cloud pressures, as well as the lack of explicit treatment for aerosols. Finally, the daily mean and median SCDs over the presumably SO2-free equatorial east Pacific have increased by only 0.0035DU and 0.003DU respectively over the entire 15-year OMI record, while the standard deviation of SCDs has grown by only 0.02DU or 10%. Such remarkable long-term stability makes the new dataset particularly suitable for detecting regional changes in SO2 pollution

Filling the Gaps: The Synergistic Application of Satellite Data for the Volcanic Ash Threat to Aviation

Although significant progress has been made in recent years, estimating volcanic ash concentration for the full extent of the airspace affected by volcanic ash remains a challenge. No single satellite, airborne or ground observing system currently exists which can sufficiently inform dispersion models to provide the degree of accuracy required to use them with a high degree of confidence for routing aircraft in and near volcanic ash. Toward this end, the detection and characterization of volcanic ash in the atmosphere may be substantially improved by integrating a wider array of observing systems and advancements in trajectory and dispersion modeling to help solve this problem. The qualitative aspect of this effort has advanced significantly in the past decade due to the increase of highly complementary observational and model data currently available. Satellite observations, especially when coupled with trajectory and dispersion models can provide a very accurate picture of the 3-dimensional location of ash clouds. The accurate estimate of the mass loading at various locations throughout the entire plume, however improving, remains elusive. This paper examines the capabilities of various satellite observation systems and postulates that model-based volcanic ash concentration maps and forecasts might be significantly improved if the various extant satellite capabilities are used together with independent, accurate mass loading data from other observing systems available to calibrate (tune) ash concentration retrievals from the satellite systems

Using Satellite Observations to Evaluate the AeroCOM Volcanic Emissions Inventory and the Dispersal of Volcanic SO2 Clouds in MERRA

Simulation of volcanic emissions in climate models requires information that describes the eruption of the emissions into the atmosphere. While the total amount of gases and aerosols released from a volcanic eruption can be readily estimated from satellite observations, information about the source parameters, like injection altitude, eruption time and duration, is often not directly known. The AeroCOM volcanic emissions inventory provides estimates of eruption source parameters and has been used to initialize volcanic emissions in reanalysis projects, like MERRA. The AeroCOM volcanic emission inventory provides an eruptions daily SO2 flux and plume top altitude, yet an eruption can be very short lived, lasting only a few hours, and emit clouds at multiple altitudes. Case studies comparing the satellite observed dispersal of volcanic SO2 clouds to simulations in MERRA have shown mixed results. Some cases show good agreement with observations Okmok (2008), while for other eruptions the observed initial SO2 mass is half of that in the simulations, Sierra Negra (2005). In other cases, the initial SO2 amount agrees with the observations but shows very different dispersal rates, Soufriere Hills (2006). In the aviation hazards community, deriving accurate source terms is crucial for monitoring and short-term forecasting (24-h) of volcanic clouds. Back trajectory methods have been developed which use satellite observations and transport models to estimate the injection altitude, eruption time, and eruption duration of observed volcanic clouds. These methods can provide eruption timing estimates on a 2-hour temporal resolution and estimate the altitude and depth of a volcanic cloud. To better understand the differences between MERRA simulations and volcanic SO2 observations, back trajectory methods are used to estimate the source term parameters for a few volcanic eruptions and compared to their corresponding entry in the AeroCOM volcanic emission inventory. The nature of these mixed results is discussed with respect to the source term estimates

Using Neural Networks to Improve the Performance of Radiative Transfer Modeling Used for Geometry Dependent Surface Lambertian-Equivalent Reflectivity Calculations

Surface Lambertian-equivalent reflectivity (LER) is important for trace gas retrievals in the direct calculation of cloud fractions and indirect calculation of the air mass factor. Current trace gas retrievals use climatological surface LER's. Surface properties that impact the bidirectional reflectance distribution function (BRDF) as well as varying satellite viewing geometry can be important for retrieval of trace gases. Geometry Dependent LER (GLER) captures these effects with its calculation of sun normalized radiances (I/F) and can be used in current LER algorithms (Vasilkov et al. 2016). Pixel by pixel radiative transfer calculations are computationally expensive for large datasets. Modern satellite missions such as the Tropospheric Monitoring Instrument (TROPOMI) produce very large datasets as they take measurements at much higher spatial and spectral resolutions. Look up table (LUT) interpolation improves the speed of radiative transfer calculations but complexity increases for non-linear functions. Neural networks perform fast calculations and can accurately predict both non-linear and linear functions with little effort

Nitrogen Dioxide Trend over the United States: the View from the Ground, the View from Space

Emissions of nitrogen oxides (NOx) are decreasing over the US due to environmental policies and technological change. We use observations of tropospheric nitrogen dioxide (NO2) columns from the Ozone Monitoring Instrument (OMI) satellite instrument and surface NO2 in-situ measurements from the air quality system (AQS) to quantify the trends, and to establish the relationship between the trends in tropospheric column and surface concentration. Both observations show substantial downward trends from 2005 to 2013, with an average reduction of 35 percent according to OMI and 38 percent according to AQS. The annual reduction rates are largest in 2005-2009: -6.2 percent per year and -7 percent per year observed by OMI and AQS, respectively. We examine various factors affecting the estimated trend in OMI NO2 columns and in-situ NO2 observations. An improved understanding of trend offers valuable insights about effectiveness of emission reduction regulations on state and federal level